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ABSTRACT 

We study the H 2 cooling emission of forming galaxies, and discuss their ob- 
servability using the future infrared facility SAFIR. Forming galaxies with mass 
< 1O 11 M emit most of their gravitational energy liberated by contraction in 
molecular hydrogen line radiation, although a large part of thermal energy at 
virialization is radiated away by the H Lya emission. For more massive objects, 
the degree of heating due to dissipation of kinetic energy is so great that the 
temperature does not drop below 10 4 K and the gravitational energy is emitted 
mainly by the Lya emission. Therefore, the total H2 luminosity attains the peak 
value of Lh 2 ~ 10 42 ergs/s for forming galaxies whose total mass M tot ~ lO n M . 
If these sources are situated at redshift z ~ 8, they can be detected by rotational 
lines of 0-0S(3) at 9.7/zm and O-OS(l) at 17//m by SAFIR. An efficient way to find 
such H2 emitters is to look at the Lya emitters, since the brightest H2 emitters 
are also luminous in the Lya emission. 

Subject headings: cosmology: theory — galaxies: formation — infrared: galaxies 

1. Introduction 

In the last decade, observational studies on primeval galaxies has made great progress. 
For example, deep surveys in the optical (Williams et al. 1996), near-infrared (Maihara et 
al. 2001) and mid-infrared (Taniguchi et al. 1997) wavelengths have revealed early evolution 
of galaxies. 
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Since galaxies in the earliest phase consist mostly of gas, knowledge on interstellar 
medium in forming galaxies is valuable for understanding galaxy formation process. Molec- 
ular line observation has been a powerful tool in diagnosing the physical conditions in local 
star- forming regions (e.g., Dyson & Williams 1997). Recently, molecular line emission has 
also been found to emanate from high-redshift galaxies and the number of detected sources 
has been increasing rapidly (e.g., Ohta et al. 1996; Omont et al. 1996; Combes 1999). 
Among theoretical studies, Spaans & Silk (1997) investigated CO line emission from proto- 
galaxies. However, the detection of metal lines depends on the degree of metal enrichment 
and thus preceeding star formation is required. 

Can we detect molecular lines in primordial pristine gas clouds? The possibility of 
detecting H 2 line emission from forming galaxies was discussed as early as the firm estab- 
lishment of the Big-Bang cosmology (Saslaw & Zipoy 1967). Other pioneering work includes 
that by Izotov & Kolesnik (1984) and Shchekinov & EnteT (1985), which treats the H 2 emis- 
sion from proto-clusters of galaxies, i.e, the so-called "Zel'dovich pancakes" . More recently, 
Ciardi & Ferrara (2001) investigated the H2 emission from the super-shells created in the 
intergalactic medium as a result of multiple SN explosions in PopIII objects. They concluded 
that about 10 per cent of the SN mechanical energy is carried away by the H 2 1-0S(1) rovi- 
brational line. Stimulated by the recent theoretical understanding of the first star-forming 
regions (e.g., Omukai & Nishi 1998; Abel, Bryan, & Norman 2002; Bromm, Coppi, & Larson 
2002; Omukai & Palla 2001, 2003), research has been vigorously carried out on the radia- 
tion from individual star-forming cores. Flower & Pineau des Forets (2001) calculated H 2 
emission from a J-type shock expected to appear in the prestellar collapse. Kamaya & Silk 
(2002) studied the H 2 emission from a collapsing prestellar core using a simplified model, 
and pointed out that the 0-0S(3) line is the strongest during the prestellar collapse phase. 
Ripamonti et al. (2002) arrived at similar conclusions by a more refined calculation. Shibai 
et al. (2001) studied the detectability of H 2 absorption lines in high-redshift objects. All 
these studies, however, suggest that the prospect for direct detection of H 2 emission from 
forming galaxies is rather dim, though not rejected, even with next generation facilities. 

On the other hand, Lya emitting galaxies in high-redshifts have been observed (e.g., 
Hu, Cowie, & McMahon 1998; Steidel et al. 2000; Kodaira et al. 2003). Stimulated by these 
findings, some authors have investigated the Lya cooling radiation from protogalaxies by 
semi-analytical modeling (Haiman, Spaans, & Quataert 2000) or by numerical simulations 
(Fardal et al. 2001). They concluded that the greater part of the gravitational energy is 
indeed emitted by the H Lya line emission. In their model, however, H 2 -line cooling is 
not included, and the temperature in protogalaxies thus remains above 10 4 K, below which 
atomic cooling is strongly suppressed. Although this is a fair assumption for the massive 
objects that they considered (see §3 below), for smaller objects with total mass < 1O U M , 
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sumcient H 2 can be formed as a result of enhanced ionization degree by virialization shock 
and its delayed non-equilibrium recombination (Shapiro & Kang 1986; Susa et al. 1998; Oh 
& Haiman 2002). By this mechanism, the temperature drops further below 10 4 K due to H 2 
line cooling. Therefore, the greater part of gravitational energy liberated by contraction will 
be available for the H 2 cooling radiation. In this paper, we study the H 2 cooling radiation 
from those forming galaxies and discuss the prospect for observing such sources. Our result 
is quite encouraging, suggesting that H 2 emission is indeed observable by the future infrared 
satellite mission SAFIR from sources as high as z = 8. 

The organization of the rest of this paper is as follows: In §2, the method of calcula- 
tion is described. In §3, the prediction for the H 2 emission from forming galaxies and its 
observability is discussed. Finally in §4, a brief summary and discussions are presented. 
Throughout the paper, we adopt a ACDM cosmology in accordance with the WMAP data 
(Spergelet al. 2003): (Qm, ^a, ^b, h, as) — (0.3, 0.7, 0.05, 0.7, 0.9). 



2. Model 

2.1. H 2 emission from forming galaxies 

In this section, we describe our model for evaluating H 2 line emission. We use a simplified 
evolution model where a forming protogalaxy is treated as a uniform cloud. 

As the initial state, we take an instance of virialization. We use the analytic formula 
for the density at virialization 

Pvir = PcritA c (1) 

where 

A c = 18^ + 82* -39,', ^- flo(1 + "* )3 + iV (2) 

which has been found to provide a good description of the numerical results (Bryan & 
Norman 1998). 

The critical density at Z — -2-vir Colli be written as 

Pcrit = Pcrit,o[^o(l + z vir) 3 + ^a], (3) 

where p C rit,o is the present-day critical density. 

The radius and temperature of objects of total mass M tot in this instance are 

\ 1/3 
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and 

T vir = §f^M tot , (5) 

respectively. 

We first describe the dissipationless case and later describe the prescription for the 
dissipation of kinetic energy. For the sake of simplicity, the dark matter density is assumed 
to be constant after the virialization. The evolution of gas density is computed by solving 
the virial equation for a spherical object 

d 2 R GM ( _ T \ 

~W\ ~f~(Kj) ( } 

where 

and M(R) is the mass enclosed within a spherical shell of radius R. 
The thermal evolution is computed by solving the energy equation 

de d / 1\ 

~V~n - - A rad - Achcm, (8) 



(9) 



dt dt \p 
where e is the specific thermal energy 

1 k K T 



7ad - 1 ' 

and A rad and A chem are the radiative and chemical cooling rates per unit mass, respectively. 

The radiative cooling rate includes H2 line emission, HD line emission, atomic emission 
by H and He, and Compton cooling. The H 2 line cooling is treated as in Omukai (2001), 
where the effect of CMB is taken into account and the ortho-para ratio of 3:1 is assumed. 
For the HD cooling, the cooling function by Galli & Palla (1998) is used. The effect of CMB 
for HD cooling is treated approximately by replacing the cooling rate by HD line Ahd(^) by 
Ahd(^) — A H D(r ra d) (e.g., Tegmark et al. 1997). The atomic emission includes collisional 
excitation cooling, recombination cooling, and bremsstrahlung of H and He. The coefficients 
given by Anninos et al. (1999) are used. 

The above formulation neglects dissipation of kinetic energy of contracting gas; after the 
initial thermal energy is carried away by radiation, essentially all the gravitational energy 
liberated by contraction contributes to accelerating the gas infall. If the gas, however, settles 
into a Kepler rotation disk, the kinetic energy is reduced by a factor of 2 compared to the free- 
falling sphere with the same radius, i.e. half of the kinetic energy must be dissipated away. 
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Such dissipation may proceed through the thermalization of kinetic energy by turbulent 
motion induced by the non-spherical nature of the collapse. To mimic this effect, we add a 
dissipation term in equations of motion and of energy; 



d 2 R GM / T 



dt 2 R 2 V T vir (R) 

and 



+ Arb (10) 



37 = -P^l ( - ) - A rad - A chcm + T trb , (11) 



dt dt \p / 

where D trh and r trb are the drag and heating terms due to dissipation of kinetic energy. 
To compute these quantities, we assume that a fraction / trb of kinetic energy is thermalized 
during the collapse. Specifically, we first evaluate at each timestep the increment of kinetic 
energy in the dissipationless case AKq using equations (6) and (8). We then solve equations 
(10) and (11) assuming that f tr bAK = r turb At is injected as heat and the actual increment 
of kinetic energy is AK = (1 — f tlh )AK . In reality, the value of / trb is likely to lie between 
0.5 (for the Kepler rotating disk) and (for the spherical collapse). We adopt / trb = 0.25 as 
our fiducial value, while we also study the dissipationless case (/ trb = 0) for comparison. 

The freeze-out values of chemical abundances are used for the initial abundances for 
the computation: y e = 3 x 10~ A ,y H2 = 1 x 10~ 6 and y HD = 1 x 10~ 9 (Galli & Palla 1998). 
The chemical reactions between H, and He compounds are solved following Omukai (2000; 
reactions 1-22), with HD-bearing reactions of D3, D4, D8,and D10 in Galli & Palla (1998) 
added. We find, however, that HD cooling is not important in all cases we studied. 

The cloud collapse is followed from the virialization until either the local Hubble time 

:i-^n) 1/2 



O-Hn 



elapses or the radius reaches the rotation barrier 



nl /2 (i + z ) 3 / 2 



(12) 



/gas \ / ^turn 
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where f gas is the mass fraction of gas. We adopt a median value 0.05 for the dimensionless 
spin parameter at turn-around A turn (Barns & Efstathiou 1987) and the gas density enhances 
by a factor of 5000 between _R vir and R rot - The duration of emission is therefore 

t cm (M, z) = mm[tu(z),t rot (M, z)\, (14) 

where t TOt (M, z) is the time taken to reach R TOt from i? vir . 
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2.2. Merger probability of protogalaxies 



In predicting the luminosity and the abundance of protogalaxies, we also take into 
account explicitly hierarchical clustering and mergers of these objects. This is of particular 
importance at high redshift, say z £ 8, where a significant fraction of galaxy-scale halos 
merges into larger halos before the free-fall time is reached. For such halos, the duration of 
emission will be shorter than that estimated in Sec 2.1. 

We incorporate the variations in the emission times by defining the "merger probability 
weighed average luminosity" of a protogalaxy with total mass M at redshift z as 



L line (M,z) 



rt(z) 
Jt(z)- 


rint ( M _ „ \dp(M,z,z f ) 
t em (M,z)^ineK IVI , Z , Z f) dz f 


dz f 
dtf 


dt f 


rt(z) dp(M,z,z f ) 
Jt(z)-t ein (M,z) dz f 


dzf 
dt f 


dt 


f 





(15) 



where Lj^ e (M, z, zf) is the intrinsic line luminosity at z of the galaxy virialized at Zf, corre- 
sponding to the cosmic time tf = t(zf), and dp(M, z, Zf)/dzf is the probability distribution 
of formation epochs of virialized halos. The former quantity is computed using the model 
described in Sec 2.1. Since t em (M,z) is typically an order of magnitude smaller than t(z), 
we neglect the changes in p v ; r and i? v ; r between z and Zf. For the latter we adopt the pre- 
scription of Lacey & Cole (1993) as described in Kitayama & Suto (1996); the formation 
epoch Zf is defined as the time at which a halo assembles at least half of its final mass. In 
the rest of this paper, unless otherwise stated, we simply refer as "luminosity" to that given 
in equation (15). 

The number density of protogalaxies in accord with the above definition of average 
luminosity is 



dn(M, z; L Une ) 



_ dn PS (M : z) f t{z) dp(M, z, z f ) 

dM dM Jt(z)-t cm {M,z) dzf 

where dn P s(M,z)/dM is the Press-Schechter (1974) mass function. 



dzi 



dt. 



dtf, 



(16) 



If strong bias is present in such high redshift we consider here, namely z ~ 10 — 20, 
the number of observable objects increases remarkably. In this sense, our result, which 
corresponds to the no bias case, can be regarded as the lower bound on the number of 
objects. 
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3. Results 

3.1. H 2 emission from individual objects 

In the following, we consider two epochs in the early universe, z =20 and 8. The former 
corresponds to the typical epoch for the first star formation (e.g., Abel et al. 2002). The 
latter is believed to be shortly before the completion of the cosmic reionization from the 
observation of high-redshift quasars (e.g., Becker et al. 2001). Direct observation of objects 
forming in this epoch is therefore intersting for understanding the reionization process. 

The total luminosity in H 2 lines is shown in Figure 1 as a function of the total mass of 
the objects, along with the Lja luminosity for the virial redshift (a)z v - ir =20 and (h)z v - w =8. 
The corresponding virial temperature and the observed flux are indicated in the top and 
right axes. The observed flux in a line is calculated using the usual relation 

Sline = L hnc /4lTd L {z) 2 , (17) 

where cIl(z) is the luminosity distance to the source at redshift z. The contribution to the 
total H 2 luminosity from each line is shown in Figure 2 for the fiducial case (/ t rb = 0.25) and 
in Figure 3 for the dissipationless case (/ tr b = 0). Important lines are listed in Table 1 (for 
z = 20) and 2 (for z — 8) for the fiducial case. 

We first focus on the fiducial case / tr b = 0.25. From the H 2 emission property, the 
protogalactic clouds can be categorized into the following four mass ranges. The temperature 
evolution for typical cases in these regimes is shown in Figure 4 for z = 20. 

. Pop III objects (M tot < 1O 8 M ) 

Small objects whose virial temperature T vir < 10 4 K cannot excite atomic hydrogen, 
and cool only by H 2 lines. Thus both the thermal energy at virialization and the grav- 
itational energy liberated by the collapse are converted to H 2 emission. The dominant 
lines are the lowest excitation ones, e.g., 0-0S(0) at 28 //m and 0-0S(l) at 17 /im (see 
Figure 2). 

• dwarf galaxies (1O 8 M < M tot < 10 n M Q ) 

Clouds with a virial temperature exceeding 10 4 K can cool by Lja radiation. As we 
can see in Figure 4, the temperature first drops rapidly from the virial temperature 
to that appropriate for H 2 cooling (~ 10 3 K). The initial thermal energy is radiated 
away by the atomic radiation processes including the Lja emission (T > 10 4 K) and by 
the H 2 line emission (T < 10 4 K). The collapse proceeds quasi-isothermally thereafter 
due to the H 2 line emission. This can be observed in Figure 1 as the emergence of 
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Lja luminosity. At the same time, the slope of the increase of the H 2 line luminosity 
toward higher mass flattens slightly around T vir ~ 10 4 K. The low excitation lines, e.g. 
O-OS(l) at 17 /iin and 0-0S(3) at 9.7 /im, are still dominant (Figure 2). 

• giant galaxies (10 n M Q < M tot < 1O 12 M ) 

Since the available gravitational energy increases toward massive objects, the heating 
rate by the dissipation of kinetic energy also increases. For objects > lO n M , the 
heating becomes so great that the temperature does not fall below 10 4 K. These objects 
collapse isothermally at 10 4 K, after initial rapid cooling (see Figure 4 a). During this 
isothermal collapse, a large fraction of energy is emitted in the Lja radiation. Then 
the Lja luminosity exceeds the H 2 luminosity (Figure 1) for those objects. Despite 
the high temperature, a small amount of H 2 is present. Only a small fraction of 
gravitational energy is emitted in the H 2 lines. This can be observed as a plateau 
in the H 2 luminosity in Figure 1, as the larger mass compensates for the lower H 2 
emissivity. Corresponding to the abrupt change in thermal property, the dominant 
lines are altered (see Figure 2). Because the collapse proceeds with high temperature, 
high excitation lines become more important. The strongest line is 0-0S(3) at 9.7/im. 
Although O-OS(l) is still in the second place, the vibrational line 1-0S(3) at 2.0/im also 
becomes important. These objects are the brightest ones in H 2 emission and the H 2 
luminosity reaches about 3 x 10 41 ergs s -1 . Note that these are also most luminous in 
Lya emission. 

• no-cooling objects (M tot > 1O 12 M ) 

For these massive objects, both Lja and H 2 luminosities decline rapidly, because, in the 
redshift range we are considering, such massive objects cannot cool within the merging 
timescale and then remain in the virialization state (Silk 1977; Rees & Ostriker 1977). 

We stress the importance of the dissipation of kinetic energy. Early studies by Saslow 
& Zipoy (1967) and by Izotov & Kolesnik (1984) did not take the dissipation into account. 
As we can see in Figure 1, the maximum H 2 luminosity would be smaller by an order 
of magnitude without the dissipation of kinematic energy. In this case, only the portion 
of the initial thermal energy below 10 4 K would be available for H 2 emission, with most 
gravitational energy being converted into kinetic energy. The discrepancy in emitted energy 
increases toward higher mass because of the larger gravitational energy to be released during 
the collapse; in the dissipationless case, the H 2 luminosity increases only linearly in mass, 
while in the dissipational case it increases approximately as oc M^. 
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3.2. Number counts of H 2 emitters and their contribution to the far-IR 

background light 

H 2 -line photons emitted at z ~ 10 have been red-shifted to the far-infrared (FIR) 
wavelength today. The future space telescope, the Single Aperture Far-Infrared Observa- 
tory 4 (SAFIR), is planned to have unprecedented sensitivity at the FIR and submillimeter 
wavelengths, with the limiting flux of ~ lCT 19 erg s _1 cm -2 in FIR. With this limiting flux, 
brightest sources at z = 8 can be observable by the 0-0S(3) 9.7/im line (Figure 2 b), while 
at redshift z = 20, even the brightest H 2 emitters are below this threshold (Figure 2 a). 

Figure 5 shows the number of sources per unit solid angle above the limiting flux of (a) 
10~ 19 erg s _1 cm~ 2 and (b) 10~ 20 erg s _1 cm -2 whose emission redshift is larger than z. For the 
limiting flux 10~ 19 erg s _1 cm~ 2 , which is the planned limit for SAFIR, although the number 
of sources with their H2 total flux exceeding this limit is about 10 6 str^ 1 as early as z = 10, 
even the strongest line 0-0S(3) 9.7 fim cannot be observed until z = 8 (see Figure 5 a). 
We can understand this from Figure 2 (b); at z — 8, although the maximum total H 2 flux 
exceeds 10 _19 erg s _1 cm~ 2 by some factor, even the most dominant line 0-0S(3) 9.7/zm barely 
reaches this limiting flux because of sharing the luminosity with other lines. 

For the limiting flux 10~ 20 erg s _1 cm -2 , the brightest objects at z = 20 could be visible 
by the 0-0S(3) 9.7/xm line and the 0-0S(l) 17/im line (Figure 2 a), although such massive 
objects are still very rare at z = 20 (Figure 5 b). For this limiting flux, the number of sources 
observable by these lines exceeds ~ 100 str" 1 at z ^ 15. 

Figure 6 further illustrates the number of H 2 emitters for each line above the limiting 
flux per unit redshift range at z — 8. At the limiting flux of SAFIR, 10~ 19 erg s _1 cm~ 2 , 
observable lines are 0-0S(3) 9.7/im, and possibly 0-0S(l) 17// m, although the latter slightly 
falls short of 10~ 19 erg s~ 1 cm~ 2 . If the limiting flux were an order of magnitude deeper, some 
vibrational lines including 1-0Q(1) 2.4 fim and 1-0S(1) 2.1/im could also be detected. 

Finally, Figure 7 displays the contribution of these H 2 emitters to the extragalactic FIR 
background radiation. Also plotted for reference is the FIR background spectrum detected 
by COBE/FIRAS (Fixsen et al. 1998). The H 2 line emission from forming galaxies makes 
only a negligible contribution to the background light, being smaller than the COBE/FIRAS 
result by several orders of magnitude. 



http://safir.jpl.nasa.gov/index.asp 
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4. Summary and Discussion 

We have evaluated the amount of H2 line radiation from forming galaxies. The H2 line 
luminosity reaches the maximum value of about 3 x 10 41 erg s" 1 for protogalaxies of total 
mass M tot ~ 1O 11 M . Below this mass scale, all the gravitational energy is emitted by H 2 
line radiation, while above this scale, a large part of the energy is emitted by Lja radiation. 
Such brightest objects at z ~ 8 are detectable by the next generation telescope SAFIR via 
the 0-0S(3) 9.7/im line emission and marginally detectable also via the 0-0S(l)17//m line. 
The cosmic FIR background radiation by the H2 emission is several orders of magnitude 
below the detected value of COBE/FIRAS, and can be safely neglected. 

The H 2 luminosity from a protogalactic cloud has also been studied by Shchekinov 
(1991). Assuming that a protogalactic cloud consists of subclumps and that the luminosity 
is given by the sum of luminosity from each subclumps, he evaluated Lr- 2 > 10 42 erg/sec from 
a protogalactic cloud with gas mass 1O 11 M corresponding to total mass M tot ~ 1O 12 M . 
On the other hand, in our model the large protogalaxies of M tot > 1O 12 M has no time to 
cool and then emit only weak H 2 radiation (total H 2 luminosty Lh 2 ~ 10 40 erg/sec). This is 
because such inhomogenous clouds are regarded as a cluster of protogalactic clouds rather 
than a single cloud in our scheme. In fact, the sum of luminosity from a cluster with total 
mass of M tot ~ 1O 12 M gives a similar value as that of Shchekinov (1991). For example, the 
luminosity of a protogalactic cloud of M tot = 1O 8 M collapsing at z = 20 is L Ha ~ 10 8 erg/sec 
(Fig. 2 a). Then the H 2 luminosity of a cluster consisting of 10 4 such clouds is ~ 10 12 erg/sec. 
Note that in order for this value to be the actual luminosity from the cluster, the collapse of 
the subclumps must be simultaneous. Since this is probably a rare occurence, this possibility 
is not considered in this paper. 

As we mentioned in §3, the brightest H 2 sources, "giant galaxies", are even more lu- 
minous in the Lja emission. This has been pointed out also by Shchekinov (1991). who 
investigated pancaking shock of primordial clouds. Consequently, the most efficient way of 
finding H 2 emitters is to look at the already detected Lya emitters, rather than looking for 
them in a blank field. The origin of Lja emitters found at high-redshift is still in dispute. If 
Lja emitters are originating in the cooling radiation, they might be bright as well as in the 
H 2 emission. On the other hand, the Lya emission might be originating in the recombina- 
tion radiation due to stellar ionization (Partridge & Peebles 1967). In this case, there is no 
reason to expect bright H 2 emission from those sources. Nonetheless, since H 2 is vulnerable 
to photodissociation by stars as we will discuss below, even if the H 2 emission will not be 
detected, this will not demonstrate the recombination origin for the Lja radiation. 

Here, we discuss some effects that can alter the amount of the H 2 luminosity from 
forming galaxies. The density contrast inside the clouds, in particular higher density than 
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average in the inner region, results in luminosity higher by some factor, reducing the col- 
lapse timescale compared to the merging timescale. On the other hand, the following effects 
reduce the H 2 luminosity: First, molecular hydrogen is vulnerable to photodissociation by 
stars. According to the recent WMAP results, the reionization epoch of the universe could 
be as early as z ~ f 7 ± 5 (Spergel et al. 2003). At least below this redshift, protogalax- 
ies are irradiated by the UV background radiation. According to our study, objects most 
luminous in H 2 emission are massive ones whose mass is about 10% of galaxies today. 
Such massive objects are not vulnerable to photodissociation by UV background radiation. 
Therefore, our results remain unchanged except for "PopIII objects" , for which the feedback 
effect is significant (Haiman, Rees, & Loeb 1997; Omukai & Nishi 1999; Ciardi et al. 2000; 
Kitayama et al. 2001; Ricotti et al. 2002). Although "giant galaxies" are not affected by the 
UV background radiation, an internal UV field created by local star formation might pho- 
todissociate molecular hydrogen and then reduce the amount of the H 2 emission. Another 
possibility for reducing the H 2 emissivity is by metal enrichment; with metals, the energy 
budget for H 2 emission is reduced by the cooling via metal fine-structure and metastable 
lines (e.g., Dalgarno & McCray 1972). For the metallicity Z > 10~ 2 Z Q , the metal line cool- 
ing dominates over the H 2 line cooling (e.g., Fall & Rees 1985). Below this metallicity, the 
H 2 cooling dominates and thus the effect of metal lines in reducing the H 2 luminosity can 
be neglected. Further study of the influence of stellar feedback and metal enrichment on the 
H 2 emission will be interesting. 

We are grateful for helpful discussions with H. Matsuo, H. Susa, and M. Umemura. 
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Figures 




Fig. 1. — The total luminosity in H2 lines and the Lya line luminosity as a function of the 
mass of protogalactic clouds that collapses at (a) z=20, and (b) z=8. The luminosity is 
averaged over the merger probability of halos as in equation (15). The corresponding virial 
temperature is also shown along the top axis. The observed flux today is shown along the 
right axis. Both the fiducial case (/ tr b = 0.25) and the dissipationless case (/ tr b = 0) are 
shown. 
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Fig. 2. — The contribution from individual lines to the H2 line luminosity in the fiducial case 
(/trb — 0.25) for the collapse redshift of (a) z=20 and (b) z=8. The thick solid line shows 
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Fig. 3. — The same as Figure 2 but for the dissipationless case (/ t rb = 0). 
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Fig. 4. — The temperature evolution of protogalaxies virializing at z = 20 as a function of 
the H number density n for (a)the fiducial case (/ tr b = 0.25), and (b)the dissipationless case 
(/trb — 0)- Illustrated are the protogalaxies of M tot = 10 7 (solid lines), 10 9 (dotted lines), and 
10 n M sun (dashed lines). 
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Fig. 5. — The number of sources below a redshift z exceeding the limiting flux (a) S\ ine > 
10~ 19 (erg s _1 cm -2 ), or (b) Sn ne > lCT 20 (erg s _1 cm~ 2 ). The thick solid line shows the total 
flux of H2 lines, and the thick dashed line, the Lya flux. The flux of individual H2 lines 
17/rni and 9.7/zm are also shown in thin lines. 
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Fig. 6. — The number count of H 2 emitting sources as a function of the minimum observed 
flux Siine per unit redshift interval at z = 8. 
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Fig. 7. — The contribution of H 2 cooling radiation to the cosmic background in the fiducial 
case (/trb = 0.25). Also shown is the IR background detected by COBE/FIRAS. 
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Table 1: The H 2 luminosity and observed 
flux from forming galaxies of total mass 
1O 7 M , 1O 9 M and 1O U M at z = 20. 
The quantities of H Ly a are also pre- 
sented. 
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Table 2: The same as Table 1 but for z = 8 



